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Approximately forty years ago it was realized that the time development of decaying systems 
might not be precisely exponential. Rolf Winter (Phys. Rev. 123, 1503 (1961)) analyzed the 
simplest nontrivial system - a particle tunneling out of a well formed by a wall and a delta-function. 
He calculated the probability current just outside the well and found irregular oscillations on a short 
time scale followed by an exponential decrease followed by more oscillations and finally by a decrease 
as a power of the time. We have reanalyzed this system, concentrating on the survival probability of 
the particle in the well rather than the probability current, and find a different short time behavior. 

PACS numbers: 03.65.-w, 03.65.Xp 



1. INTRODUCTION 



Unstable quantum mechanical systems are generally described in terms of an exponential decay law: the probability 
P(t) that a given system remains after a time i, called the survival probability, is given by 

P(t) = e- Xt , (1) 

where A is a constant. However it has been known since Khalfin's work Q in 1958 that decays could not be pre- 
cisely exponential. At very long times P(t) must decrease more slowly than Eq. ([j]) and, given somewhat restrictive 
assumptions on the energy distribution of the system, the derivative of P(t) must vanish at t = 

In 1961 Winter gave a concrete example of nonexponential decay by analyzing a simple one-dimensional tunneling 
problem. He calculated the probability current just outside a delta-function barrier to determine the rate of decay. 
Winter gave approximate analytic expressions for the wave function, but to be precise he was forced to calculate the 
current numerically. Given the state of computers in 1961 even to calculate the current was a computational tour de 
force and he did not calculate the survival probability itself. 

Until recently, the direct experimental observation of departures from exponential decay has proved elusive. Indeed, 
precise (3 decay experiments covering a range of decay times less than lO -4 ^^ in 60 Co and out to 45^/2 in 56 Mn 
showed no deviation from exponential decay [|J. More recently, several experiments have reported measurement- 
induced suppression in quantum systems at short times. One || observed the Zeno effect ||, which is a result of 
dP(t)/dt\t=o being zero, in stimulated emission. Another Q examined the effects of interaction-free measurements 
with photons. In the past few years, the first experiments on an unstable quantum system jsj, |9| were performed. 
These experiments measure the time dependent Landau-Zener tunneling of trapped sodium atoms. They observe 
nonexponential decay directly and also see a Zeno effect and an anti-Zeno effect. 

In view of the fact that some observations are now being made, we have revisited Winter's calculation. We believe 
it is basically correct. Here, we simply extend it slightly by calculating the survival probability which is the quantity 
of direct experimental interest. Some of the nonexponential effects in the probability current, such as large oscillations 
at short times, are considerably damped in the survival probability. Other effects, such as large oscillations during 
the transition from exponential to power law fall off at large times, are not damped. We turn now to the details. 



2. WINTER'S MODEL 

In this model, the potential consists of an infinite potential barrier at x = —a and a delta-function at x = 0, 

I oo for x < —a 
I Uo (x) for x > —a 
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with U > 0. The initial state is taken to be totally within the well formed by the barrier and the delta-function 

(2/a) 1 / 2 sin(n7rx/a) for — a < x < 

for x < —a or x > 



^(x,0) 



(3) 



To find the time dependence of -0, we solve Schrodinger's equation for the energy eigenstates <f>k{x). Taking into 
account the vanishing of the wave function at x — —a, we can write 



(t>k(x) 



A sinfc(x + a) for —a < x < 
Be lkx + Ce~ ikx for x > 



(4) 



where k = (2mE) 1 ^ 2 with m and E denoting the mass and energy of the particle in the well. The boundary conditions 
at x — then relate the constants A, B and G as 
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for —a < a; < 



fx) = A 



where G = 2maU and q — ka. This then gives 

sin(q(x/a + 1)) 

sm(q(x/a + 1)) + ^ sin(q) sin(qa;/a) for a; > 
with A determined by normalization. We choose to normalize the (f>k (x) as 

/oo 
dx <j)* k i (x)(j)k{x) = 6(fc' — £;) . 
-a 

In view of Eq. (Q) , we have 

dx(j)* k ,{x)4>k{x) = 



1 H sin q cos q H sin q cos q H sin q sin q 

9 9 99' 
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ax cos k a; sin kx -\ — - 



q Jo 
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dx cos kx sin k s 



dx sin k! x sin fca; 



sinq' sinq > \ A\ 



Using the integrals 



dx sin k'x sin kx — — fc) — <5(fc' + fc)] 



1 f°° If 00 

— / da; cos k'x sin fca; + — / dx cos kx sm k' x = ir 2 S(k)5(k') , 
k Jo k J 

and the fact that k, k' > 0, Eq. ([^) gives 



\A\ 



Ti" sjq 2 + qG sin 2q + G 2 sin 2 q 
Expanding t/>(a;, 0), Eq. (||), in terms of the (j>k{x) using Eq. jl^), the expansion coefficients C(k) are 

nnA _ 2ny/wa qsinq 

r - ^9 + qG sin 2q + G 2 sin 2 q 

The wave function ip(x,t) for i > can then be expanded as 



ip(x,t) = / dkC{k)(/) k (x)e- 



2n 



dq- 



-iTq 2 - 



qsinq /(a;,q) 



(5) 
(6) 

(7) 

(8) 



q 2 - n 2 TT 2 q 2 + Gq sin 2q + G 2 sin q 



(9) 

(10) 
(11) 

(12) 
(13) 

(14) 
(15) 
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where T = t/2ma 2 . In terms of t = x/a, the function f(x,q) is 



f gsin(£+l)<7 for — a < x < 

f( X ' q ^ = \ ■ to , ,\ , n ■ -Of ^ n ' ( 16 ) 

I q sm(t + l)q + G sin g sinfg tor x > 



3. DISCUSSION 



Using the wave function of Eq. (|lq), Winter calculated the probability current at the barrier, j(x,T) at x = 0, 

(17) 



ax ax 



We repeated this calculation and our results, shown in Fig.[l], agree with his. In particular they show large oscillations 
at T < 0.8. Like Winter, we used G = 6 and n = 1. 

However the survival probability 

P(T)= f dx\^(x,T)\ 2 (18) 



has only small oscillations, on the order of 5% as shown in Fig.|| and Fig.[| This is not too surprising since, by the 
continuity equation, P(T) and j(0,T) are related by 



P{T) - P(0) = -2ma 2 [ dT 1 j(0, T' 
Jo 



) (19) 

and the effect of the integral, in this case, is to smooth out the oscillations. 

The source of the oscillations is the second denominator of Eq. ([l5]), which gives an infinite number of poles in the 
complex g-plane. Closing the contour in such a way as to pick up the contributions from these poles gives a sum of 
exponentials with different energies and lifetimes 

e -^ 2 ^ e -i Si T e -T/2n _ (20) 

Winter gives approximate analytic expressions for the pole positions as a function of G but these expressions are not 
accurate for G < 20. In Table| we show the pole positions for several values of G as well as the lifetime and energy 
from each pole. Since the first lifetime is significantly larger, the other terms die out quickly but not before causing 
the deviations shown in Figs. || and ||. In fact the difference of the first two energies gives a period of oscillation of 
0.25 which agrees exactly with these figures. 

Once these small short time oscillations die out, the slope determined numerically and shown in Figs.|| and |] agrees 
precisely with that given by the first pole in Table ffl. This shows that there are no additional contributions from the 
extra pieces of the contour required to close it until we get to large times. At large times the factor exp(— iTq 2 ) 
oscillates rapidly and the contribution to the integral is negligible except for q very small. As a consequence, we can 
expand the integrand in powers of q and the integral in Eq. ( |l5| ) goes as 

^" iT9 V ~ ^ (2D 

More precisely, if we close the contour at an angle of —45° with the real axis the extra contribution, beyond the poles 
discussed above, is 

* { ^' 2 n Aotm^Tm + °{M' - a - xi "' (22) 

So at very large times the survival probability goes as 1/T 3 . This agrees with our numerical results as shown in Fig.^ 
for T greater than 18. 

As the time dependence changes from exponential to an inverse power Winter found large oscillations in the 
probability current. Indeed those also occur in the survival probability as shown in Fig.||. They arise from the 
interference between the time behavior shown in Eq. (|2l]) or Eq. (^||) , which have no oscillatory dependence, and that 
of Eq. (20). The energy from Table|, e = 7.59, gives a period of oscillation of 0.82 which agrees exactly with the 
numerical results of Fig.|i[ 
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4. SUMMARY 



For this simple example we can understand the decay through a barrier at all times. At small time the survival 
probability has small oscillations because of interference between exponentials with different lifetimes. This is followed 
by an exponential region once all the exponentials except the one with the longest lifetime decay away. Oscillations 
occur during the transition to the 1 /T 3 behavior which holds at large time. All this is summarized in Fig. ^| Except 
for the reduced magnitude of the small time oscillations these results for the survival probability are essentially what 
Winter found for the probability current. 

We close with a few comments. 

The size and duration of the small time oscillations are limited in this system by the largest lifetime being signifi- 
cantly longer than the others. However one could conceive of a situation where the longest two or more lifetimes are 
similar in magnitude. Then the oscillations could last longer and be larger, perhaps even to the extent that there 
would be no single pure exponential region. While it is easy enough to construct a C(k) which will give this behavior, 
it is not clear that such a function is associated with a physically realizable system. 

We define the survival probability as the probability the particle remains in the well. Usually the survival probability 
is defined in terms of the survival of the initial state 

P(T) = |W0)|V(T))| 2 (23) 

2 

dxip*(x,0)ip(x,T) . (24) 

-a 

Since Winter used the probability current just outside the barrier it seems clear our definition agrees with his. For 
this problem it makes very little difference which definition is used. This is illustrated at short times for G = 6 in 
Fig. ^J, where the definition Eq. (p4|), computed using the first two poles in Table [i], is compared with the corresponding 
two pole evaluation of Eq. (|lj). The two pole approximation turns out to be quite accurate for < T < 2. We have 
checked that this result holds for other values of G. For long times, where we can use Eq. (p2|), the ratio of Eq. (|l8| ) 
to Eq. ( p4] ) is 7r 2 /6. The combination of the two pole approximation for the short time behavior of the survival 
probabilities Eqs. (|l^) and (^) together with the long time behavior given by Eq. ((2^) is shown in Fig.^j. 

Note that the wavefunction, Eq. (Q), must be equal to the initial wavefunction, Eq. (||), at t = and therefore the 
current at the barrier, x — 0, is zero at t = 0. Alternately, simply substitute the explicit form for ip( x i t) into Eq. ( |l7j ) 
and note that each surviving term goes as sin(q 2 T). Thus 

|p(0) = -i(0,0)=0 (25) 

in agreement with more general arguments 

Since the small time behavior of P{T) is larger than pure exponential (see Fig.|2[) the Zeno effect holds. The survival 
probability at a particular time, say T — 1.0, will be larger if intermediate measurements are made. For example, if we 
compare ten measurements at intervals of T = 0.1 and two measurements at intervals of T = 0.5 with no additional 
measurements then, using P(T = 0.1) = 0.9119, P(T = 0.5) = 0.4754 and P{T = 1.0) = 0.2194, we find 

[P(T = 0.1)] 10 > [P(T = 0.5)] 2 > P(T = 1.0) . (26) 

At very large times where P(T) is again not exponential there is an anti-Zeno effect - intermediate measurements 
make the survival probability much smaller || . For example 

[P(T = 15)] 2 < P(T = 30) . (27) 

There can also be a short time anti-Zeno effect where the decay rate is enhanced by repeated measurements at a 
time interval chosen to occur near the minimum of the oscillations |l(], [ll], 12, ljj]. Observation of this is reported in 



Such an effect is system dependent and may not occur at all. Indeed, for our canonical case of G = 6 it does not 
occur (note that the in Fig. || are always greater than unity). However, we have checked that it does occur for other 
values of G, specifically G = 20. 

Finally it should be noted that there are deeper questions concerning quantum mechanics in a time-asymmetric 
universe that we have not considered here fll4|. 



5 



Acknowledgments 

One of us (D.A.D.) thanks Arno Bohm and Xerxes Tata for teaching him the basics of nonexponcntial decay. This 
research was supported in part by the National Science Foundation under grant PHY-0070443 and by the United 
States Department of Energy under Contract No.DE-FG03-93ER40757. 



[1] L.A. Khalfin, Soviet Phys. JETP 6, 1053 (1958); Soviet Phys. Doklady 5, 515 (1960). 

[2] J.J. Sakurai, Modern Quantum Mechanics (Addison- Wesley, New York, 1994). Also see L. Fonda, G.C. Ghirardi, and A. 

Rimini, Rep. Prog. Phys. 41, 587 (1978). 
[3] R.G. Winter, Phys. Rev. 123, 1503 (1961). 

[4] E. B. Norman, S. B. Gazes, S. G. Crane and D. A. Bennet, Phys. Rev. Lett. 60, 2246 (1988). 
[5] W. M. Rano et. al, Phys. Rev. A41, 2295 (1990). 

[6] B. Misra and E.C.G. Sudarshan, J. Math. Phys. 18, 756 (1977); C.B. Chiu, E.C.G. Sudarshan and B. Misra, Phys. Rev. 
D16, 520 (1977). 

[7] P. Kwiat, H. Weinfurter, T. Herzog, A. Zeilinger and M. A. Kasevich, Phys. Rev. Lett. 74, 4763 (1995). 

[8] SR. Wilkinson et. al, Nature 387, 575 (1997); Q. Niu and M. G. Raizcn, Phys. Rev. Lett. 80, 3491 (1998). 

[9] M. C. Fischer, G. Gutierrez-Medina and M. G. Raizen, Phys. Rev. Lett. 87, 040402 (2001). 
[10] W. C. Schieve, L. P. Horwitz, and J. Levitan, Phys. Lett. A 136, 264 (1989). 
[11] A. G. Kofman and G. Kurizki, Phys. Rev. A 54, R3750 (1996). 
[12] A. G. Kofman and G. Kurizki, Nature (London) 405, 546 (2000). 
[13] P. Facchi, H. Nakazato, and S. Pascazio, Phys. Rev. Lett. 86, 2699 (2001). 
[14] See A. Bohm, Phys. Rev. A60, 861 (1999) and references therein. 



Re(g) 


Im(g) 


r 


e 


G=l 


2.2986 


-0.76605 


0.142 


4.70 


5.4340 


-1.19699 


0.038 


28.10 


8.5857 


-1.42451 


0.020 


71.69 


11.7351 


-1.57997 


0.013 


135.2 


G=6 


2.7579 


-0.14043 


0.646 


7.59 


5.7135 


-0.37015 


0.118 


32.51 


8.7753 


-0.55535 


0.051 


76.70 


11.8767 


-0.69721 


0.030 


140.6 


G=20 


2.9958 


-0.02054 


4.063 


8.97 


6.0109 


-0.07438 


0.559 


36.13 


9.0532 


-0.14565 


0.190 


81.94 


12.1204 


-0.22148 


0.093 


146.9 



TABLE I: Positions of the poles in the integrand of the expression for ip(x,i). For each value of G the positions of the first 
4 poles are given, as well as the lifetime, r = — [4Re(g) Im(g)] -1 , and energy, e = (Re(g)) 2 — (Im(g)) 2 , for each pole. In the 
figures we use G = 6. 
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FIG. 1: The deviation of the probability current Eq. ([n]) from exponential behavior at short times using r = 0.646. 
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FIG. 2: The deviation of the survival probability Eq. from exponential behavior at short times using r : 



0.646. 
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Survival Probability 
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FIG. 3: The survival probability at short times. 
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FIG. 4: The survival probability at all times. 
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FIG. 5: The short time survival probability calculated from Eq. ( |24[ ) (solid line) is compared with the survival probability calculated from 
Eq. (Il8h (dashed line) in the two pole approximation. 




